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ABSTRACT 

The accuracy requirements on station-keeping for constellations of free-flying collectors cou- 
pled as (future) imaging arrays in space for astrophysics applications are examined. The basic 
imaging element of these arrays is the two-element interferometer. Accurate knowledge of two 
quantities is required: the projected baseline length, which is the distance between the two in- 
terferometer elements projected on the plane tranverse to the line of sight to the target; and 
the optical path difference, which is the difference in the distances from that transverse plane to 
the beam combiner. "Rules-of-thumb" are determined for the typical accuracy required on these 
parameters. The requirement on the projected baseline length is a knowledge requirement and 
depends on the angular size of the targets of interest; it is generally at a level of half a meter 
for typical stellar targets, decreasing to perhaps a few centimeters only for the widest attainable 
fields of view. The requirement on the optical path difference is a control requirement and is 
much tighter, depending on the bandwidth of the signal; it is at a level of half a wavelength 
for narrow (few %) signal bands, decreasing to « 0.2 A for the broadest bandwidths expected 
to be useful. Translation of these requirements into engineering requirements on station-keeping 
accuracy depends on the specific details of the collector constellation geometry. Several examples 
are provided to guide future application of the criteria presented here. Some implications for the 
design of such collector constellations and for the methods used to transform the information 
acquired into images are discussed. 



Subject headings: 

1. Introduction 

One of the major problems affecting the de- 
sign of future systems for high-resolution astro- 
nomical imaging using constellations of free-flying 
collectors in spac^ll is the necessity to precisely 
maintain the "optical figure" or surface accuracy 
of the equivalent aperture for extended periods 
of time. On sufficiently-bright targets, this accu- 
racy might be achieved using signal photons from 



Some examples of such systems presently under study by 
NASA include the "Stellar Imager" (SI), the "Terrestrial 
Planet Finder - Interferometer" (TPF-I), and the SPECS 
mission study for a sub-millimeter space interferometer. 



the target itself in order to operate various con- 
trol loops. Such "adaptive" control systems can 
be found, for instance, driving deformable mirrors 
on filled-aperture telescopes, and controlling de- 
lay lines on interferometers. However, in the more 
general (and often more interesting) case where 
the target is very faint, the accuracy requirement 
translates into a tight requirement on the geome- 
try of the optical system. In the case of a constella- 
tion of free- flying collectors in space, the accuracy 
requirement on the optical figure of the instrument 
for observations of faint targets becomes a require- 
ment on station-keeping. 

Station-keeping may be either inertial, with re- 
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spect to some fixed reference points (such as dis- 
tant quasars), or relative to the rest of the space- 
craft in the constellation. Relative station-keeping 
may be done so as to yield a good optical fig- 
ure, but still leave the whole constellation tum- 
bling about some arbitrary axis, so this alone is 
not sufficient. A significant design effort has been 
expended on dealing with the "inertial" part of the 
problem, and most approaches for measuring the 
overall rotation of the whole constellation involve 
the use of the " rotat ing Sagnac interferometer" 
described e.g. byl Hecht (2002 jl. 

Early studies l|Jonesl Il99ll : IESAI Il996l ) of the 
station-keeping problem with free-flying collectors 
took as a "straw-man" criterion that the positional 
accuracy required would be of the order of a small 
fraction of a wavelength, e.g. w A/50. This trans- 
lates to « 10 nm at an optical wavelength of 500 
nm. A combination of radio and laser-ranging 
systems may be adequate for station-keeping at 
a level of microns, but more elaborate measures 
will be required in order to reach the nanometer 
level if such accuracy is really required. Current 
thoughts for achieving this level of accuracy in the 
absence of sufficient photons from the target itself 
include fine control using photons from another 
nearby bright "phase reference" star in the field of 
view, or the development of major technical im- 
provements in ranging systems. 

The use of a phase reference star has been de- 
veloped extensively in order to overcome severe 
fluctuations in the fringe pattern caused by at- 
mospheric turbulence in ground-based optical and 
near-IR long- baseline inter fe rome ters (see e.g. the 
summary by iQuirrenbach ( 2000l )). Application 
to the (more liesurely but analogous) problem of 
station-keeping in a slowly-wandering space inter- 
ferometer constellation was proposed e.g. in an 
ESA study more than ten years ago ( ESA 19961 ) 
where it was shown that the tolerances in the di- 
rections perpendicular and parallel to the direc- 
tion to the target could be reduced by factors of 
order 10'^ and 10®, respectively, if a bright phase- 
reference star existed within an angular distance 
of 2' from the target of interest. However, the 
penalty here is a reduction in that fraction of the 
celestial sphere which can be observed, since suit- 



able phase-reference stars may not be available for 
every target of interest. In addition, the effective 
field of view of the constellation's optical system is 
likely to be very smalU, thereby further reducing 
the utility of this approach. 

Compared to the ground-based atmospheric- 
turbulence problem, the longer time scales ex- 
pected for station-keeping drifts may make it feasi- 
ble to use photons from the target itself to stabilize 
the fringes. This can be done in certain special- 
ized applications (e.g. nulling interferometers such 
as TPF-I and Darwin), but becomes increasingly 
difficult in general for faint targets. In addition, 
the target should not be appreciably resolved at 
the interferometer baselines in use, since the sig- 
nal photon level then also decreases. But the main 
science goal may well be to measure just this kind 
of structure in the target; hence the need to use 
target photons for fringe stabilization is in general 
very restrictive. 

Since we can expect that current laser-ranging 
technology will continue to improv^, one can ask 
whether it might be possible to relax the « A/50 
requirement enough to accomplish station-keeping 
purely by on-board laser-ranging methods. This 
would permit the observation of arbitrarily-faint 
targets located anywhere in the sky. Questions 
which arise include: How firm is the "A/50" re- 
quirement? Does the physics of producing im- 
ages from interferometer fringes justify this num- 
ber? Are there ways of extracting the desired im- 
age information, perhaps only approximately, but 
which can tolerate larger errors? The search for 
answers to these questions has motivated the work 
described in this paper. 

An example of the kind of space-based imaging 
system considered here is shown in Figure [1] This 
is a cross-section sketch of a concept design for the 
Stellar Imager, a constellation of about 30 free- 
flying collectors operating at optical/UV wave- 
lengths and distributed over an area with maxi- 
mum dimension of up to 500 meter^. These col- 
lectors direct the radiation they receive to a cen- 



^One design using this concept is called the "Kilometric Op- 
tical Gyro". 



•^For example, this would be the case with certain designs 

for the Stellar Imager (Figure [T}. 
^Impressive gains in precision laser metrology have recently 

been achieved for use in the NASA/JPL Space Interferom- 

etry Mission (SIM). 
^The Stella r Imager mission study is described at: 

|http: / /hires. gsfc.nasa.gov/ ~si| 
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Fig. 1. — Side view sketch of a possible design for the "Stellar Imager" concept (not to scale). The collectors, 
which are identical 1-meter spherical mirrors in this design, reflect the light from the target (off to the far 
right in the figure in the direction indicated by the arrow "L") to a central hub spacecraft. 



tral "hub" spacecraft for further processing and 
downlink. Figure [5] shows a possible distribution 
of the collectors as seen from the central hub. 

The plan of this paper is as follows. First, a 
simple (but entirely general) model of the imaging 
process for free- flying constellations of collectors is 
described. The basic imaging element of these con- 
stellations is the two-element interferometer, and 
the image quality is particularly sensitive to errors 
in the determination (and stability) of the phase 
of the interferometer fringe pattern. This model 
leads to the identification of two important pa- 
rameters which play a role in determining the ac- 
curacy with which the imaging can be done. The 
first of these is the projected baseline length, which 
is the distance between the two collector elements 
of each constituent interferometer projected on a 
plane perpendicular to the line of sight to the tar- 
get, and the second is the optical path difference, 
which is the difference in the distances from that 
perpendicular plane through each collector to the 
interferometer beam combiner. 

At this point it is useful to clarify the usage 
here of the terms knowledge of a parameter and 



its control. In general, precise control of baseline 
length is not required; what is required is pre- 
cise knowledge of the actual values at any time. 
That knowledge permits us to place a data sample 
at the right location in the aperture plane before 
carrying out the imaging computations. On the 
other hand, precise knowledge of the optical path 
difference is not sufficient; we must control this 
quantity to be close to zero at all times. Just how 
this knowledge and control will be provided will 
depend on the design of the particular constella- 
tion; those details actually do not concern us here, 
but they will affect the distribution of the error 
budget over the various parts of the instrument. 
What does concern us here is to determine the 
required accuracy on knowledge of the projected 
baseline length and control of optical path differ- 
ence in order to achieve a given level of accuracy in 
the measurements which will go into the final syn- 
thesized image. As we shall see, the physics of the 
interferometric imaging process implies that the 
typical knowledge accuracy required for the pro- 
jected baseline length Bfj turns out to depend on 
the angular size of the target of interest; it is gen- 
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"Stellar Imager" Concept Mission 
View of collector constellation from direction of target 




500 m 



Fig. 2. — View of the "Stellar Imager" collector constellation from the target showing one of the many 
interferometer baselines which make up this telescope "aperture" . This view is in the transverse plane 
(indicated by the arrow "T" in Figure [1]), and is the distance between any two collectors projected on 
this plane. 



erally at a level of half a meter for typical stellar 
targets, decreasing to a few centimeters only for 
the widest practical fields of view. On the other 
hand, the control accuracy required for the opti- 
cal path difference OPD is much higher, and de- 
pends on the bandwidth of the signal. It is at a 
level of half a wavelength for narrow (few %) signal 
bands, decreasing to « 0.2A only for the broadest 
bandwidths expected to be useful. This is a fac- 
tor 10 less stringent than the rough value of A/50 
used in the initial exploratory studies mentioned 
earlier, but exploiting this relaxed value will re- 
quire post-processing of the interferometer data in 
a computer. 

The paper continues with a brief discussion 
on translating the two parameters identified into 
more familiar terms such as the range and bearing 
from one spacecraft to another. This discussion 
is sketchy and given only by way of illustration. 



since the details of specific constellations are not 
the main topic here. Finally, some implications 
of these results for the future design of free-fiyer 
imaging systems are discussed. Conventional "di- 
rect" imaging systems, which form images directly 
on a panoramic detector (e.g. a CCD) at the focus 
of a ("dilute") aperture, do not easily permit the 
use of the relaxed requirements on station-keeping 
described here. "Indirect" imaging systems, where 
the images are formed by post-processing of inter- 
ferometer data in a computer, provide opportuni- 
ties to make use of this (and other ancillary) infor- 
mation, and it is likely that future high-resolution 
space-based imaging systems for astronomy will 
be of the indirect type. 

Background material on the topics of this paper 
can be found in several excellent tutorials on the 
application of interferometry techniques to optical 
astronomy. For example, see the volume edited 
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by iLawson (20o3il. Among others, reviews of 
pro gress and recent advanc es have b een pubhshed 
by Shao k, Colavital ( 1992 ) and by iQuirrenbach 
()200lh . 



2. The Imaging Process 

How do we make images using data collected 
from astronomical targets with a constellation of 
free-flying collectors? The relevant physics is de- 
scribed mathematically in terms of the mutual 
coherence function of the wavefront, a quantity 
which can be measured with an interferometer, 
and the Van Cittert-Zernike theorem, which re- 
lates this mutual coherence function to the ob- 
ject on the sky (i.e. the target brightness distri- 
bution) by a Fourier transform. The re are numer- 
ous te x tbooks on th is topic, see e.g. iBorn &: WoH 
(|l975l ): iHechtl (|2002l ). The result relevant for now 
is that, whether the imaging system is "direct", 
where an image is instantaneously formed on a 
2-dimensional detector (e.g. a dilute-aperture im- 
ager), or "indirect", where an image is formed 
by post-processing in a computer (e.g. an array 
of Michelson interferometers), the image is built 
up by summing fringe patterns formed by the in- 
terference of light gathered from all pairs {i,j) of 
collector elements in the constellation. One such 
pair, connected as a Michelson interferometer, is 
sketched in Figure [H 

At any given wavelength A, the interferome- 
ter fri nge pattern has the form (e.g. iBorn &: WoH 
()l975l) . Chapter 7): 



= Po[l + Acos(27rfcA)], 



(1) 



where Pq is the total signal (photons. Joules, etc.) 
collected by the two elements, k — 1/X, A is the 
fringe amplitude, (0 < ^ < 1), and A is the optical 
path difference (often designated as OPD) com- 
posed of: 



Bij sin ( 



(2) 



A is the difference in total path (meters) from the 
wavefront plane through the two "arms" to the 
beam combiner, consisting of Bij smO + di on the 



^Lawson also maintains the Optical Long Baseline Interfer- 
ometry News web site with a wealth of useful links and 
reference information at; 

|http://olbin.jpl.nasa.gov/papers/index.litml| 



left side of Figure [3l and dj on the right. The 
angle 9 (in radians, often assumed small so that 
sin 9 ~ 9) is measured in the plane containing the 
interferometer baselin^ll- As we shall see momen- 
tarily when we consider detection systems with fi- 
nite bandwidth, it is important to keep A « 0, 



so if 



we shall need to have di 



This 



can be achieved by bouncing the light around, in 
space or in a delay line, before it enters the beam 
combiner. For instance, if the constellation is de- 
ployed on a parabolic surface and the beam com- 
biner is located at the focus of this parabola, then 
the delays for each collector are "automatically" 
adjusted to be (nearly) equal. This particular ge- 
ometry is called a "dilute aperture" since it re- 
sembles a conventional reflecting optical telescope 
with most of the collecting area removed. But 
this particular geometry is a special case and is 
not a requirement; the only requirement is that 
the optical paths somehow be kept nearly equal. 
We also define the fringe period measured on the 
sky as X/Bfj radians. Note that the fringe pe- 
riod increases as the projected baseline decreases 
with increasing 9. Note also that a typical system 
at optical wavelengths may have Bf^/X « 1 x 10^ 
and 6* w 3 - 100 x IQ-i" rads, or « 60 - 2000 ^as. 

It is convenient to simplify the way we write 
the fringe pattern by defining the quantities (fix = 
kBijSin9 and cjjj = k(dj — di), both in units of 
"turns", so that Equation [1] becomes: 



P^J =Po{l + ^COs[2^(0T "</-/)]}, 



(3) 



where (f>T is called the fringe phase and depends 
on the location of the target on the sky, and (j>j is 
the instrument phase, which can be controlled by 
altering the internal delay difference dj — di . The 
elemental narrow-band fringe pattern of Equation 
[3] is sketched as a function of instrument phase in 
the top panel of Figured] Note that the structure 
of the fringes can be determined by scanning the 
instrument phase (f)i over one or two turns and 
fitting the resultant fringe pattern to Equation 
[3l thereby obtaining the total power, the fringe 
amplitude, and the fringe phase. This form of 



^This direction is in fact on the surface of a cone with half- 
opening angle of (7r/2) — 6 and axis coincident with the 
baseline Bij. The "ideal" interferometer has a constant 
response on the surface of this cone, and the field of view 
must be further restricted e.g. by adding an aperture stop 
into the optical system. 
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Fig. 3. — The basic geometry of an interferometer. Two collectors i and j sample the wavefront from a very 
distant target at two points separated by Bf^ = Bij cos 6 meters. The two samples are correlated in the beam 
combiner after experiencing relative path delays of Bij sin 9 + di and dj , respectively. The beam combiner 



calculates the fringe parameters, e.g. by scanning S — dj 
fitting the pattern. 



di over one or two fringe periods, sampling, and 



the equation also makes it clear that errors in 
the instrument phase </>/, such as those arising 
from station-keeping errors, directly masquerade 
as changes to the measured value of fringe phase 
(/)T, and therefore to errors in the measured posi- 
tion 9 of the target. Finally, Figure |4] also makes 
it clear that the position of the target in the field 
of view is obtained by counting fringes (and frac- 
tions thereof) from a reference point, for example 
the direction where 0/ = 0. 

The discussion so far has assumed that the in- 
terferometer is sensitive to only one wavelength, 
i.e. that the bandwidth of the detection system 
is infinitely narrow. The fringe pattern for a fi- 
nite bandwidth detection system is the weighted 
(vector) sum of the fringe patterns at all con- 
stituent wavelengths, with amplitudes specified by 
the product of the source spectrum and the instru- 
ment response at each wavelengtf(f|. A simple but 



illustrative case which is easy to calculate is for a 
rectangular bandpass of central wavelength Aq and 
width SX, and a target spectrum that is flat over 
this wavelength range. In this case the resulting 
fringe pattern is: 



Poll 



sin(7rA/Lc) 



X Acos(27rfcoA)], (4) 



'This expression can be written as a Fourier Transform. 



where the symbols not already defined in the pre- 
vious two equations are kg = I/Aq, and Lc = 
Xq/6X. This latter quantity is called the coher- 
ence length for t he light beam, a con cept which 
arises as follows (|Born fc Wolj (|l975l ). Ch. VII; 
Hecht (2002), Ch. 7): For a finite detector band- 
width of Sv hertz, the signal in free space is a 
compact packet of many interfering waves with a 
coherence time of r = 1/(5 seconds. Express- 
ing the bandwidth in terms of wavelength, we 
have Sv = — (c/Aq) x 6X, so the coherence time 
T = (— 1/c) X (Ag/(5A), and cr is a scale size of 
the wave packet called the coherence length, Lc = 
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Fig. 4. — Top panel: Fringe pattern from Equation [3] f^or baseline i? as a function of fringe phase 0t 
in "turns" (see text). The instrument phase is assumed to be at ~ 0. Bottom panel: Fringe pattern for 
a basehne which is in reahty 5% larger. Notice that the error grows with angular distance from the point 
where (/>/ = 0, reaching w 0.15 fringes (54° « 1 rad) at the edge of this small 3-fringe field. The y-axis is in 
units of the signal strength obtained by a single collector. 



CT — Ag/JA. The lower panel of Figure [5] shows 
a sketch of this situation for <5A/Ao — 0.2. Note 
that, within a few fringes of the central peak, the 
finite-bandwidth fringe in the lower panel shows 
the same periodicity (and positioning, or fringe 
phase) as the single- wavelength case in the upper 
panel. This means that the relative placement of 
different components of the reconstructed target 
brightness distribution will be correct, but their 
amplitudes may be wrong. Perhaps even more se- 
rious, the signal-to-noise ratio of the amplitude 
measurement will be worse. 

We are now in a position to begin the discussion 
on the accuracy with which the relative positions 
of the collectors need to be known. For complete- 
ness, note that we are here focussing only on that 
part of the error budget which arises because of 
errors in the station-keeping. A full discussion of 
all error sources and their various contributions is, 
of course, dependent on the design details of each 
specific constellation. 



3. Two Requirements 

Inspection of Equations[l][2l and[4]above makes 
it clear that, in order to apply these equations to 
real data, there are two parameters which depend 
on station-keeping and which are important for 
each interferometer in the constellation: the length 
of the projected baseline Bfj, and the optical path 
difference A at the beam combiner. The accuracy 
requirements on these two parameters are pre- 
sented below first by stating the result, then pro- 
viding the derivation, and finally discussing some 
of the implications. 

3.1. Baseline 

The fractional accuracy SBfj/Bfj with which 
the projected distance between any i,j collector pair 
in the transverse plane must be known is: 

SB/B <l/{NfQT). (5) 

where Nf is the number of fringes in the angle p, p 
is the angular radius of the field of view containing 
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the object being imaged, and Qt is the "transverse 
quahty factor" expressing the accuracy X/Qt with 
which the "optical figure" of the imaging system 
is to be held. The super- and subscripts of B have 
been dropped for convenience. 

3.1.1. Demonstration: 

The fringe pattern for an elemental interferome- 
ter expands and contracts about the field center as 
B decreases or increases owing to station-keeping 
errors, and for a given (but in reality unknown) 
baseline error of SB the resulting fringe phase er- 
ror grows with angular distance from the center of 
the field. The situation for the narrow-band case 
is sketched in Figure [H 

Without loss of generality we can take the cen- 
ter of the fringe pattern to be the centroid of 
the target brightness distribution at the operat- 
ing wavelength A. Let p be the angular radius (in 
radians) of a circle on the sky which is just large 
enough to include all the emission from the target. 
The number Nf of fringes in p is Nf = p/(\/B). 
This number changes with small changes in the 
baseline length according to 5Nf = p5B/\. Sub- 
stituting for p from the initial expression for Nf 
results in 6Nf = Nf x SB/B. If we require that 
the optical quality of the signal provided by this 
interferometer be better than a fraction 1/Qt of 
a wavelength, then we must have 6Nf < 1/Qt- 
From this we derive SB/B < \/ {N f Qt), which 
proves the result. 

3.1.2. Discussion: 

As an example, suppose the source is a distant 
star and the instrument is intended to provide 30 
independent resolution elements across the stellar 
diameter. Each resolution element is of angular 
size \/Bmax radians, where Bmax is the longest 
(projected) basehne in the constellation. Suppose 
further that it can be safely assumed there is no 
signal outside a larger circle of some given angu- 
lar radius, say p = 20 x X/B,nax radians on the 
sky measured from the center of the star. For this 
example, Nf w 20. If we demand an "optical qual- 
ity" of, say, A/50 in our imaging instrument, then 
Qt = 50 and dB/B < 1/1000. For a maximum 
baseline of i? = 500 meters, the largest unknown 
SB we can tolerate is « 0.5 m. Note that since 
Nf (X B for a given source size and operating 



wavelength, SB is actually independent of B, so 
the required accuracy in this example is the same 
for every baseline in the constellation. However, 
the stringency of the requirement grows with the 
size of the target field of view in the sky. For 
instance, for a supergiant star the desired FOV 
could be a factor of 10 larger than the example 
calculated above, and the required tolerance on 
any baseline would be « 5 cm. But now such a 
constellation would provide 300 independent res- 
olution elements across the star. If the science 
goals permit, it may be more prudent to "shrink" 
the overall extent of the constellation by a factor 
10 in order to return to the scale of the initial ex- 
ample and relax the station-keeping requirements 
accordingly. 

It is interesting to ask what the limit on SB / B 
might be for very large fields of view. For instance, 
for a Michelson radio interferometer with a sin- 
gle "feed horn" at the focus of each collector, the 
largest field of view will be limited by the diffrac- 
tion pattern of the individual collector elements. 
This field has an angular radius « A/a where a 
is the diameter of a collector element. In this ex- 
treme case, Nf « i?/a, so the maximum tolerable 
error is: 

SB/B<l/{NfQT)^a/{BQT); (6) 

i.e., SB < u/Qt, which depends only on the size 
of the elements in the constellation. If a typical 
collector is one meter in diameter and Qt = 50, 
the required knowledge accuracy is w 2 cm. The 
situation gets more demanding if a wide field 
beam combiner is used, such as an optical Michel- 
son interferometer with detector pixels that are 
larger than the collector diffraction pattern, or 
with many independent pixels in the detection sys- 
tenl!|. 

To conclude this section, note that changes in 
operating wavelength have similar effects on the 
array reponse to those discussed above and can be 
modelled in the same way. One particularly dele- 
terious effect of this is that the point spread func- 
tion of a system with a finite spectral resolution 
will vary over the field of view, becoming radially 
elongated at large angular distances. Such effects 



Dilute aperture "Fizeau" designs are not as optimum in 
signal-to-noise as are the Michelson designs, and will not 
be discussed here any further. 
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were first analysed in detail by radio astronomers 
modeling the reponse of grou nd-based r adio syn- 
thesis imaging arrays (see e.g. I Thompson. (,1994 )). 

3.2. Optical Path Difference 

The accuracy 6A with which the optical path 
difference must be controlled to zero is: 

6A/X0 < {1/Qopd) X (Xo/SX). (7) 

where SX/Xq is the fractional bandwidth of the 
signal and Qo pd is the "OPD quality factor" , by 
analogy to the Qt factor defined earlier. This re- 
quirement can also be written as SA < Lc/Qopd, 
where Lc is the coherence length defined following 
Equation [H 

3.2.1. Demonstration: 

Some clarification is needed first, since the real 
problem here may not be entirely obvious. It is 
clear from Equation [1] that errors of 5A/X directly 
translate to a fringe error of G 5A/X where G is a 
geometry-dependent factor of order unity depend- 
ing on the specific geometry of the constellatiorF^I. 
However, in this case it is not an expansion or con- 
traction of the fringe pattern about the field cen- 
ter, but a gross translation of the entire fringe pat- 
tern in a direction parallel to that of the baseline 
projection on the field of view. At first sight this 
is a serious problem, since it appears to require 
holding the collector positions to a tiny fraction 
of the wavelength. However, since it is the whole 
fringe pattern that moved and not the relative sep- 
aration between two points in the field of view, the 
angular distance measured in fringes between the 
center and the edge of the field remains the same. 
In other words, the relative fringe phase is actu- 
ally not sensitive to this shift! Figure [5] shows a 
sketch of the situation. 

If we can design our system and/or (post-) 
process the data in order to use this "relative" 
phase insead of an "absolute" phase, we would be 
entirely insensitive to OPD errors, although we 
would lose the ability to measure the target po- 
sition in an absolute sense. Except for accurate 
wide-angle astrometry, no other applications in as- 
trophysics require such absolute position informa- 
tion. The real problem comes about because of a 



For example, G ~ 2 for a classical (Fizeau) imaging system. 



loss of fringe amplitude at large non-zero values of 
the OPD owing to the finite bandwidth of the de- 
tection system. Only if the bandwidth is infinitely 
narrow will the fringe pattern amplitude remain 
independent of the OPD error 5 A as shown e.g. 
in the upper panel of Figure [6l the lower panel of 
this figure shows a "real" fringe pattern for a finite 
(20%) bandwidth following Equation HI It is clear 
that, as the OPD error grows, so does the error 
in the measurement of the fringe amplitude. The 
scale size of the fringe pattern envelope is char- 
acterized by the coherence length Lc as shown in 
Equation [H If we can keep the OPD error 5 A to 
some fraction 1/Qopd of the coherence length Lc, 
then this amplitude error can be minimized and 
the measurement of fringe amplitude and fringe 
phase can be made subject only to photon noise 
in the time scale of the drift in the station position. 
This requirement becomes 5 A < {1/Qopd) x Lc 
or 5A/Xo < {^/Qopd) x (Ao/(5A), which proves 
the result. 

3.2.2. Discussion: 

As an example, for Qopd = 50 and a band- 
width of a few %, say Xq/SX « 30, the control 
requirement on SA is 0.6 A. This requirement will 
still be difficult to achieve at optical wavelengths, 
but now there is a better chance that on-board 
laser ranging systems may eventually provide ad- 
equate position knowledge without the use of sig- 
nal photons, thereby allowing long integrations on 
faint sources. Bandwidths in excess of « 10% in 
a single detector channel are not likely to be used 
since the spectral dependence of the fringe pattern 
is often an important diagnostic for the on-board 
fringe detection algorithms. This means that the 
worst case for the control accuracy of the OPD 
(which is essentially the station-keeping require- 
ment, as we shall see below) is w 0.2 Aq, still a 
factor of 10 less stringent than the early estimates 
of 0.02 Aq mentioned in the introduction to this 
paper. 

Now that the two major physical quantities 
have been identified, and their accuracy require- 
ments estimated, we can consider how to apply 
these criteria to real cases of specific mission con- 
cepts. 
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Fig. 5. — Top panel: Narrow-band fringe pattern from Figured) Bottom panel: Narrow-band fringe pattern 
shifted by a phase error of S(j) — I rad, caused e.g. by an OPD error of SA « A/6. 



4. Applications to Specific Mission Con- 
cepts 

The previous section has identified the two 
physical parameters B and A which characterize 
the constituent fringe patterns that go into build- 
ing an image from data obtained with a constella- 
tion of free- flying collectors. In order to take the 
next step of translating the general knowledge re- 
quirements on SB and SA into nanometers and mi- 
croradians for a specific mission concept, we need 
to distinguish two cases depending on the bright- 
ness of the targets to be observed. 

4.1. Bright Targets 

If enough photons will be available to operate 
a servo loop to keep the optical path difference 
A « to within some fraction of the coherence 
length Lc, then the control requirement on SA 
need not translate into a station-keeping require- 
ment, but rather into a requirement on the total 
stroke and accuracy of the on-board delay line. In 
this case we are left only with the knowledge re- 
quirement on the projected baseline length which 
depends on the size of the target field of view. 



as stated in Equation [51 The requirement on the 
accuracy in the range from one spacecraft to an- 
other is of order half a meter. This still permits 
imaging of the stellar surface with a relative reso- 
lution of about 30. An error in the bearing of the 
second spacecraft with respect to the first in the 
interferometer pair also appears as an error in the 
projected baseline, but this component will be less 
significant; from the geometry of the interferom- 
eter in Figure [3] it is easy to show that an error 
of Sa in the bearing translates into a fractional 
projected baseline error of SB/B « (fo)^/2. For 
the example taken in gUl 6B/B < 1/1000, so the 
requirement on the bearing is 0.1 rad « 1° inde- 
pendent of B . 

4.2. Faint Targets 

If we want to observe faint (or partially re- 
solved) targets which do not provide enough pho- 
tons to properly operate a control loop on the op- 
tical path difference, then Equation [7] drives the 
station- keeping requirement. Errors in station- 
keeping map into errors in A in a manner depen- 
dent on the specific constellation geometry. For 
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Fig. 6. — Top panel: Fringe pattern as a function of OPD error (in wavelengths) for a very narrow-band 
signal. Bottom panel: The same fringe pattern for a square bandpass of 20%. Note how the fringes 
disappear beyond Nf fa 5. 



instance, for a design where the beam combiner 
lies between the two collectors (as in the TPF-I 
concept), a station-keeping error of SB in the base- 
line length will produce a delay error of 6 A = 6B, 
and we have seen from the demonstration follow- 
ing Equation [7] that such errors need to be kept 
smaller than some fraction 1/Qopd of the coher- 
ence length. If laser ranging is available and capa- 
ble of measuring the baseline length to sufficient 
precision, then that information could be fed to 
an on-board delay line in the target signal path to 
correct for baseline length errors. This would turn 
the requirement on SB from a control requirement 
into a knowledge rcg uiremenlF^. 

Bearing errors 6a also directly translate into 
delay errors; for the SI concept of Figure [1] this 
is SA « BSa, so the knowledge of bearing of 
the second spacecraft with respect to the first is 
5a < Lc/ (BQopd) which for sample values used 
in the discussion in H3.2\ becomes Sa < 0.6 X/B 
rads or « 0.1 milli-arcseconds for B = 500 m and 



^^However, the required precision on that knowledge remains 
a fraction 1/QoPD of th^ target signal coherence length. 



A — 500 nm. This will be a difficult requirement 
to meet if the collectors in the constellation are 
all more-or-less in a plane. Placing an additional 
"metrology" spacecraft some distance from this 
plane (e.g. in the "Hub" spacecraft for the SI de- 
sign of Figure [1]) would alleviate the requirement 
by converting the inter-collector bearing measure- 
ment for station-keeping errors in the direction to 
the target (the "L" direction in Figure [T]) to a 
range measurement (with the usual required ac- 
curacy of order Lc/Qopp) from a collector to the 
metrology spacecrafl^. 

5. Indirect Imaging 

In order to take advantage of the relaxed con- 
trol requirement on the optical path difference de- 
scribed in Equation [71 we must find a method 
of turning the recorded fringe amplitudes and 
(slowly-drifting) phases into an image in such a 



Note that the problem remains of rendering the station- 
keeping inertial, in the sense described in the Introduction 
to this paper. The presence of a "hub" spacecraft may 
present an opportunity to address that problem as well. 



11 



way as to be insensitive to shifts of the whole fringe 
pattern owing to small OPD errors. There are 
several related "post-processing" techniques which 
could be used for this. For instance, if the field of 
view contains a bright star besides the target of in- 
terest, fringe phases could be "referenced" to this 
star. Another approach (first used at radio wave- 
lengths) is to reference phases to a bright spectral 
line feature. 

One quite general approach to this "unstable- 
phase" imaging problem which has been used with 
success for « 40 years in the radio astronomy 
community uses the c oncept of "closu re phase". 
Almost 50 years ago, Jennison ( IQSSh presented 
a technique for measuring relative fringe phase 
which used three radio-linked collectors coupled as 
three interferometers operating at a wavelength of 
2.4 meters over baselines up to « 10 km. Owing 
to a variety of instrumental problems related to 
the amplifiers and local oscillator electronics avail- 
able at the time, the fringe phase between any 
two collectors was unstable and could normally 
not be measured; the source structure information 
had to be derived from the (squared) fringe am- 
plitudes alone (a familiar situation in present-day 
ground-based optical interferometry). Jennison 
showed that, if the three observed fringe phases 
were summed, the resultant combined phase was 
insensitiv e to equipment i nstabilities. With this 
approach. I Jennison fc Lath am (1959) showed that 
the brightness distribution of the radio source 
Cygnus A, which until then was only known to 
be elongated, actually consisted of two separated 
sources of nearly equal brightness straddling a pe- 
culiar optical object tentatively identified at the 
time with two galaxies in collision. This was the 
first observation to reveal the double-lobed struc- 
ture of a powerful radio galaxy. 

Applications of this method to circumvent at- 
mospheric phase instabi lities in optical in terferom- 
etry were described by Ijennison' (1961*) and, ap- 
parently independently, by Rogstad (1968i) . The 
first use of the term "closure phase " seems to be 
in the paper bv lRogers et.al.l ( 1974 ) describing an 
application at radio wavelengths using very stable 
and accurate, but independent, reference oscilla- 
tors at the three stations in a so-called "very-long- 
baseline" interferometer array. Since that time, 
closure phase has been used extensively at radio, 
IR, and optical wavelengths, and there are many 



papers descri bing the subjec t, its virtues, and its 
limitation onmen (|200 3h gives a review of op- 

tical interferometry in astronomy which includes 
a discussion of closure phase, and also points out 
the intimate relationship between it and the com- 
plex bispectrum, a quantity constructed in process- 
ing optical speckle interferometry dat J^. For the 
type of collector constellations discussed in this 
paper, the phase of the complex bispectrum is pre- 
cisely the closure phase, and by (vector) averag- 
ing the complex bispectrum the closure phases can 
be obtained in principle for arbitrarily faint tar- 
gets. One important point to make here is that 
use of such techniques requires that the fringe sig- 
nals from each pair of collectors be separable one 
from another, such that the individual interferom- 
eter fringe phases can be uniquely assigned to the 
appropriate pair of coUectorfJij. This leads in turn 
to requirements on the design of the optics at the 
location where the collector beams are brought to- 
gether to interfere. 

6. Conclusions 

The requirements on station-keeping for con- 
stellations of free-flying collectors coupled as (fu- 
ture) imaging arrays in space for astrophysics ap- 
plications have been discussed. The typical knowl- 
edge accuracy required on the projected baseline 
of the instrument depends on the angular size of 
the targets of interest; it is generally at a level of 
half a meter for typical stellar targets, becoming of 
order a few centimeters only for the widest attain- 
able fields of view. The typical control accuracy 
required on the optical path difference depends on 
the wavelength and the bandwidth of the signal, 
and is at a level of half a wavelength for narrow 
(few %) signal bands, becoming « 0.2 A for the 
broadest bandwidths expected to be useful. If the 
fringe detection system is able to use signal pho- 
tons or photons from a nearby reference star, then 
the most stringent requirement on station-keeping 



I'^For newcomers, the lectures presented at the Michelson 
Summer School by John Monnier and David Buscher in 
2001, and by Peter Tuthill in 2003 are good sources; see 
http://msc.caltech.edu/school/2003 and links there. 

^ ^Cornwelll l|l987l ) specifically (and elegantly) addresses the 
relationship between speckle masking and phase closure. 

^^Note there is no requirement here that the set of baselines 
available in the constellation be "non-redundant" , although 
that may be useful for other reasons. 
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is that on the baselines, as described in If ob- 
servations on faint targets are required, and fringe 
tracking is "blind" , then the most stringent re- 
quirement on station-keeping is that on the optical 
path difference described in H3.'2[ If on-board laser 
metrology is available and can provide knowledge 
of the internal OPD to a fraction of the target sig- 
nal coherence length, then active fringe tracking 
could compensate for station-keeping errors, even 
for faint sources. 

In any case, the requirements are less severe 
than assumed in the early studies of the problem. 
The significance of this result is that, at this re- 
laxed level of accuracy, it may be possible to pro- 
vide the necessary knowledge of array geometry 
without the use of signal photons, thereby allow- 
ing observations of arbitrarily-faint targets. 

Such constellations of free-flyers will produce 
images using various computer-based image recon- 
struction techniques which relax the level of ac- 
curacy required on the fringe phase stability of 
each component interferometer. "Closure-phase" 
imaging is one such technique which has been 
very successfully applied in ground-based radio 
and optical astronomy in order to surmount in- 
stabilities owing to equipment and to the atmo- 
sphere. This technique appears to be directly ap- 
plicable to space imaging arrays, where station- 
keeping drifts play the same role as atmospheric 
and equipment instabilities. More detailed model- 
ing of anticipated station-keeping errors and their 
effects on synthetic images is needed. For instance, 
there is at present little justification for the values 
"Qopu — Qt ~ 50" used repeatedly in this pa- 
per. A study needs to be done of the degradation 
in image quality as these parameters are dimin- 
ished (and the OPD and baseline errors grow), as 
well as a comparison of the merits of various im- 
age reconstruction algorithms including the use of 
closure phase and bispectrum analysis. 
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